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In this article we study the electrical conductivity in QED using the resummed 3PI effective 
action. We work to 3-loop order in the effective action. We show that the resulting expression 
for the conductivity is explicitly gauge invariant, and that the integral equations that resum the 
£ — \ pinching and colinear contributions are produced naturally by the formalism. All leading order 

f*"^ ■ terms are included, without the need for any kind of power counting arguments. 
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IT) ! I. INTRODUCTION 

It is well known that selective resummations can play an important role in quantum field theory. A prominant 
Qh example is the hard thermal loop theory, which includes screening effects that regulate infra-red divergences. The 
convergence of such perturbation theories can typically be improved using an n-particle irreducible (nPI) effective 
action. In addition, nPI approximation schemes are of particular interest because they can be used to study far from 
, equilibrium systems. 

The nPI effective action is consistent with the global symmetries of the theory, however, the ward identities may 
not be satisfied during intermediate steps of the calculation [2, [2] ■ To address this problem we look at the resummed 
■ nPI effective action, which is defined with respect to the self consistent solutions of the n-point functions. This 
\Q ' type of strategy was originally proposed by Baym and Kadanoff [j| and has been discussed in the context of scalar 
theories in Ref. The renormalization of resummed theories has been discussed in Ref. The resummed 

effective action respects all symmetry properties of the theory, and the n-point functions which are obtained by 
functional differentiation satisfy the ward identities. In this paper we study the resummed qed 3PI effective action. 
0^ ' We demonstrate how the calculation of transport coefficients is organized in the framework of this formalism. 

'Pinch singularities' play an important role in the calculation of transport coefficients. The basic idea is that there 
is an infinite number of terms that all contribute at the same order because of the low frequency limit in the kubo 
formula (Eqn. (J])). This limit produces pairs of retarded and advanced propagators which carry the same momenta. 
When integrating a term of the form J dpo G ret (P)G adv (P) , the integration contour is 'pinched' between poles on 
each side of the real axis, and the integral contains a divergence called a 'pinch singularity.' These divergences are 
regulated by using resummed propagators which account for the finite width of thermal excitations. This procedure 
introduces extra factors of the coupling in the denominators which change the power counting. As a consequence, 
there is an infinite set of graphs which contain products of pinching pairs that all need to be resummed. This set of 
graphs is resummed by solving an integral equation whose kernel is the square of the matrix elements that correspond 
to the 2 — > 2 scattering and production processes. 

For the qed electrical conductivity, it has already been demonstrated that this integral equation can be obtained 
from the 2PI formalism. In Ref. it was shown that the 2-loop contribution to the 2PI effective action produces 
the square of the s-channel, which gives the complete result at the leading-log order of accuracy. In a previous paper 
we have shown that the 3-loop terms in the 2PI effective action contribute the missing t- and u- channels so that 
the full matrix element corresponding to all binary scattering and production processes is obtained. 

However, the result obtained in Ref. [13] is not complete at leading order. In gauge theories, in addition to pinch 
singularities, the presence of collinear singularities makes the 1 — ► 2 scatterings as important as the 2^2 scatterings. 
These collinear terms must also be resummed by solving an integral equation. The complete leading order result is 
then obtained by solving these two coupled integral equations. In this paper we derive both of these integral equations 
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directly from the 3-loop resummed 3PI effective action. Thus we show that the full leading order contribution to the 
qed electrical conductivity can be obtained directly from the 3PI formalism. 

Our result agrees with that obtained previously in [13T |. This calculation is not obtained directly from quantum 
field theory but is derived from kinetic theory. In addition, the conductivity has recently been calculated using 
a diagramatic method based on power counting arguments in which the ward identity is used explicitly to select 
contributions that will produce a gauge invariant result Using our method, gauge invarience is automatically 

satisfied, and all leading order terms appear without the need for any kind of power counting arguments. The method 
we develop in this paper should be generalizable to the calculation of other transport coefficients, like the shear 
viscosity. 

Our calculation provides a field theoretic connection to the kinetic theory results of [13], which is useful in itself. In 
addition, it seems likely that quantum field theory provides a better framework than kinetic theory for calculations 
beyond leading order. Our results provide strong support for the use of nPI effective theories as a method to study 
the equilibration of quantum fields. 

This paper is organized as follows. In section[n]we define some notation. In section [Til] we review the 3PI formalism. 
In section [IV] we discuss the resummed 3PI formalism. In IIV Al we give the equations of motion, in IIVBI we derive 
the corresponding bethe-salpater equations, and in IIV CI we obtain an expression for the external photon propagator 
that satisfies the usual ward identity. In section [V] we use these results to obtain an expression for the conductivity 
that is complete to leading order. In section [VI] we present our conclusions and discuss future directions. 



II. NOTATION 



The electrical conductivity can be obtained from the kubo formula: 

a = U^-2Imp ii (q ,0)) (1) 
6 \dq J qo ^o 

p ll (x,y) = (iW(y)> ; fix) = ^(xW(x) 

We can write the conductivity in terms of the polarization tensor using: 

P ii (?o,o) = -4n;U(?o,o) (2) 

e 

In the rest of this paper, the coupling constant e will be absorbed into the definition of the vertex functions, and not 
written explicitly. 

In order to simplify the notation we will use a single numerical subscript to represent all continuous and discrete 
indices. For example: a photon field is written A\ :— (x); the fermion propagator is written S12 := S^Jxx, X2), etc. 
We also use an einstein convention in which a repeated index implies a sum over discrete variables and an integration 
over space-time variables. We define 2-point green functions and n-point vertex functions as shown in Fig. [1] 




FIG. 1: Definitions of notation for propagators and vertices. 



III. THE 3PI FORMALISM 



The 3PI effective action to 3-loop order can be written (l5| : 

T[tl>,$,A,S,D,V,U] (3) 
= SelfoklM] + ^TrLn^ 1 + ^Tr [(D^y 1 (D 21 - D° 21 )] - iTrLnS^ - *Tr [(S° 12 y 1 (S 21 - S° x )] 

+ r° [s, d,v,u} + +r 2 nt [s, d,v,u] 

where S is the fermion propagator, D is the photon propagator, V is the fermion-photon vertex, U is the 3-photon 
vertex (which is zero at the tree level), S c i[ip, -0, A] is the classical action, and So and Dq are the free propagators 
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given by: 



s 2 s c 



(^r 1 = 



S 2 S C 



6ip 2 6ipi SA 2 5A X 
The functions T%[S, D, V, U] and Tf^S, D, V, U] are shown graphically in Fig. H 



r, u = i 





•i/2 • 4 +i/12 4 • +i/3 




+i/4 




-i/24 



(4) 



FIG. 2: 3PI effective action to 3-loop order. 

The first diagram in Fig. [2] contains a bare fermion-photon vertex which is defined as 

S 3 S cl 



V132 — i 



51p 2 5Asstp\ 



(5) 



The equations of motion of the mean fields, propagators and vertices are obtained from the stationarity of the 
action: 



ST[^,A,S,D,V,U] =Q Sr[il;,$,A,S,D,V,U] = Q . Sr[iP,j,, A, S, D,V,U] =Q 



SA 

Sr[^,A, S,D,V,U] 

ss 

5Y[i>,j),A, S,D,V,U] 
SV 



Sip 

5rty,$,A,S,D,V,U] 
SD 

ST[^,A,S,D,V,U] 
5U 



Si' 



(6) 



IV. THE 3PI RESUMMED EFFECTIVE ACTION 

The last four equations in ((6]) can be solved simultaneously for the self-consistent solutions: S[tjj,^j,A], D[ijj , ■ip , A] , 
V[ip,ip,A} and U[ip,tp,A] which satisfy: 



STfafaA, S,D,V,U] 



ss 

5T[tl>,ji,A,S,D,V,U] 



5T[il>,$,A, S,D,V,U] 



S, D, V, U 



sv 



SD 

6r[ip,$,A,S,D,V,U] 



su 



s, d, v, u 



s, d, v, u 



S, D, V, U 

Substituting the self consistent solutions we obtain the resummed action: 

f[V>, i>, A] = T[i, i, A, Sty, i>, A],D[tP, i>, A], V[tp, A], Uty, & A}] 
The equivalence of ((3j) and ® at the exact level was shown in [16]. 
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A. Constraint Equations 

Now we look at the equations of motion ([7]). The first two equations have the form of dyson equations: 

.ST 2 [S,D,V,U] 



$12 — 0^12) 



'12 



£12 = -i- 



D12 = (D^y 1 - n 12 ; n 12 = 2% 



SS21 

ST 2 [S,D,V,U] 



SD21 



S, D, V, U 



S, D, V, U 



(9) 
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where we have defined T 2 = T® + T^ 11 *. We substitute J3j into ([9]) and take the derivatives. The results are given in 
Eqn. (9). For the fermion-photon vertex the order of indices is always (outgoing fermion, photon, incoming fermion). 
In each term, the factors that carry fermion indices are written in square brackets, and the order of the indices reflects 
the flow of fermion number. The order of the indices in the photon propagator and the 3-photon vertex does not 
matter, because of the symmetry of these n-point functions. The results are: 

§12 = (^r 1 - HViwS Vi Vh 2 ] ■ D 33 , + liVi^SviV^] ■ D 33 , (10) 

+ [Vi5'e' Se'eVe 3 'T St7Vy54' S4'4V4 3 2] ■ D^/D 33 i + \Vi4'T StyV^^'q' Sq'qVq 3 2] • U 3 i4§D 33 i D441 D55/ 
Dl2 = (^12) 1 + ^[^4:13' S 3 ' 3 V 324 , S4'4\ — i[V413' S 3 ' 3 V 3 24' S4'4] + U\ 3 > 4U 234' D 3 , 3 D 4* 4 

— [V514'S , 4'4V47'3'S'3'3V326'<S'6'6V675''S'5'5] ' £>77' — 2[Vsi4' S4'4V43'7' SV'7V765'<S'5'5] • U2 3 & D 33 > Dqq' 
+ 2^14'5^3M7'^5'67^236'^3'3^4'4-D5'5-C ) 6'6£'7'7 

These equations are shown in Fig. [3j 




FIG. 3: Self energy and polarization tensor. 
Using (f3j the second two equations in J7]) can be written: 

Vl32 = V& 2 + iVl7vh>^SwV &T 2D V7 + iV7 5ll D 5 , a U 536 D 66 ,V26>rST7 (11) 
^132 = —2iVYi 5 >S 5 > 5 V 53 QS 6e /V e ' 27' St 7 + it/715/ D 5 i 5 U 53e D 6e ,U e i27' &T7 

These equations can be represented graphically as shown in Fig. 21 
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A 



FIG. 4: Self-consistent equations for the vertices. 



B. Bethe-Salpeter Equations 



Taking derivatives of the dyson equations we obtain a set of bethe-salpeter type integral equations. In order to 
write these equations in a simple way, we define the following vertex functions: 



A - 



^r 1 

SA 3 



A132 — 



5A 3 



O - 1 Sf) 12 



(12) 



These functions are always written with the external photon on the middle leg, and this external photon is distinguished 
by an arrow in figures. To clarify the notation, all of the types of vertices that we use in this paper are shown in Fig. 
1 



= -iV 



- =-iA 



= -iU 



» - =-iQ. 



FIG. 5: The vertices. 



Some additional useful relations can be obtained from the identities: 

#13 <§32 = S12 ; D^D 32 = S12 (13) 
Differentiating (fT3|) with respect to A and using lfT2|) gives: 

^ll = S lv A V32 ,S 2 , 2 ; S ^- = 2D 11 ,fl V32 ,D 2 , 2 (14) 
6A 3 0A 3 

The bethe-salpeter equations are obtained by differentiating Q and using the definitions lfT2|) to obtain: 

A - A° + SYl12 • O - 1 61112 hK) 
A 132 -A 132 + — , n 132 _-_ (15) 

where £12 and II12 are obtained from Eqns ® and (fT0]l (and shown in Fig. [3]). When taking derivatives, numerous 
cancellations occur between contributions from different pieces of the self energy. To show how these cancellations 
occur, we distinguish two different kinds of propagators: 'rung propagators' do not connect, at either end, to a vertex 
that attaches to an external leg, and non-rung propagators do connect to a vertex that attaches to an external leg. 
For example, the third diagram in Fig. [3] gives a contribution to £ that has one fermion rung propagator, two fermion 
non-rung propagators and two photon non-rung propagators. After differentiating and collecting terms, the result is 
that all 2-loop terms containing derivatives of non-rung propagators cancel. We illustrate this point by looking at 
contributions to 62^12/ SA 3 containing derivatives of non-rung photon propagators. We give the results in terms of 
diagrams only. Differentiating the first two contributions to £ in Fig. [3] and using fTI)! gives the terms shown in Fig. 

El 



FIG. 6: Result from differentiating non-rung photon propagators in the first two contributions to — i£ in Fig. [3] 
We rewrite this result by substituting in for Vq using Eqn. l[TTj) (shown in Fig. [4]) and obtain the terms shown in Fig. 

m 



-4 



FIG. 7: Result after substituting in for Vo- 
Differentiating the third and fourth contributions to — iT, in Fig. [3] we obtain the terms shown in Fig. [HI 



FIG. 8: Result from differentiating non-rung photon propagators in the second two contributions to — iS in Fig. 

Combining the results shown in Figs. [7] and [U we see that all 2-loop terms containing derivatives of non-rung photon 
propagators cancel. The same effect occurs for derivatives of non-rung fermion propagators. In Fig. [9] we give the 
result for the first equation in (fT5| to 2-loop order and in Fig. [TOlwe give the result for the second equation in (fT5| to 
1-loop order. In section [V] we will show that higher loop terms do not contribute to the conductivity at leading order. 




FIG. 9: Bethe-salpeter equantion for the A vertex to 2-loop order. 



I 
I 

+ 




V 

I 
I 

t 



FIG. 10: Bethe-salpeter equantion for the f2 vertex to 1-loop order. 
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C. External Propagator 



The external photon propagator is defined as 



SA 2 5A 



-f[^,A] 



(16) 



In this section we show that the propagator defined in (fl6| satisfies the usual ward identities. The basic mechanism 
is simple: the dyson equations ((9]) contain s-channel ladder resummations and the bethe-salpeter equations lfl5|) 
introduce t- and u-channels, and thus restore the crossing symmetry. 

We start by taking the derivative of the modified effective action using the chain rule. In the rest of this section 
we suppress the arguments and write T[ip, tp, A, S, D, V, U] as T. We use the notation Xi to indicate one of the set 
X := {S, D, V, U} and X t to indicate one of the set X := {S, D, V, U}. The sum \ Yl&j gives the 6 terms in the set 
of pairs Xi ^ Xj where order is not respected. We write the result: 



(Dlf) 



6 2 r 



SA 2 5A 1 



x 



E 



5 2 T I 5Xi sx t 



s 2 t 



6X 2 \xSA 2 5A l 

5Xi i ^ s 2 r 



(17) 



SXiSAi \xSA 2 



E 



SXiSXq 



SXi 5Xj 
x 5Ai 5A 2 



{1-2} 



This expression can be further simplified by using the set of equations obtained by differentiating the equations of 
motion: 



6 


\ 6T 




SA 


I SXi 


X- 







5 2 T 



x 



E 



s 2 r 



SXjSXi 



SXj 
x~5A 







SXiSA 

We solve these four equations for the four derivatives: S 2 T /8Xf\y and substitute into (fl7| to obtain: 

6 2 T lfv^ S 2 T 



(18) 



SA 2 SAi 

Expanding the sum and using (f!4|) we have: 

<5 2 r i r 8 2 r 



x 



He 



SXiSAx 



SXj 
xSA 2 



{1-2} 



(Dtf)- 1 



Using |(3]) we obtain: 



6A 2 SA! 

s 2 r 



SS^SAi 



5U 3i5 5A 1 



SU: 



345 



X 

6 2 r 



(5 , 33'^-3'24''S'4'4) 
SV 3A5 



8 2 T 



x 8A 2 SV 3i5 5 A x 



x SA 2 



SDuSAi 
{1-2} 



x 



(-D33'^3'24'-D4'4) 



(19) 



(20) 



S 2 T 



SA 2 SAi 

5 2 T 
SS^SAj. 

6 2 T 



(21) 



,. ^3° 4 )- 1 
' 5A X 
5 2 T 



iA 







314 

S 2 T 



= 



which gives the final result: 



6D 3i 5Ai 6V345SA1 SUusSAi 



(Dtf)- 1 = (D^)- 1 + zA° 14 • {S M ,A 4 , 23 ,S y3 } 



(22) 



This expression has the form of the usual schwinger-dyson equation. Following the method of [12] it is straightforward 
to show that the external photon propagator in Il22|) satisfies the ward identity. 
From Eqns. (J9j) and i!22l) we extract the vertex part of the 2-point function: 



TTCXt 

11 12 



* (^314^44' A4> 23 >S 3 i 3 ) 



(23) 
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This result is illustrated in Fig. ITTl 




FIG. 11: Graphical representation of Eqn lf23)l . 



V. LEADING ORDER CONDUCTIVITY 

The conductivity is obtained from the polarization tensor through Eqns. |T]) and |(2|). The polarization tensor 
depends on the propagator S and the vertex A through Eqn. ([23| . The two equations in (fTO)) . the two equations in 
(fTTjl . and the equations represented in Figs. [UandQIil form a set of six coupled integral equations for the six quantities 
{5, D, V, U, A, fi}. In this section we show that these equations contain all of the leading order contributions to 
the conductivity. 

The pinch singularities are resummed by an integral equation of the form: 

A132 =A° 32 + J2 iM iiU ■ [&5' A 5'34' S4/4] (24) 

jE{a,b,c,d,e,f,g} 

This equation is shown schematically in Fig. [321 

y~ ,= y • -■> ••••• 

FIG. 12: Graphical representation of Eqn lf24]l . 

(i) 

The functions M^.^ represent 4-point functions where the variables on each side of the semicolon indicate legs that 
join pinching pairs of propagators. We show below that the seven terms included in the sum over j £ {a, b, c, d, e, /, g} 
correspond to the seven diagrams in Fig. [TH We start from the equation for A shown in Fig. [9] and substitute 
iteratively the equations for D, S, U, V and fl, collecting all 2-loop contributions. We describe this process step by 
step. 

Substituting the equation for 17 (Fig. [101) m t° the equation for A (Fig. [9]) and keeping terms up to 2-loop order, 
we obtain the result shown in Fig. [T3l In Fig [14] we redraw the figure in a more conventional form. 




FIG. 13: Equation for SI substituted into the equation for A. 




FIG. 14: Fig. [T3l redrawn. 
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The third and fourth diagrams on the right hand side of (fT4| are already 2-loop order. The fourth diagram is shown 
in part (c) of Fig. [151 We rewrite the third diagram as the sum of the two terms shown in parts (f) and (g) of Fig. 
fl5l Note that in Fig. [15] we include some arrows to indicate the direction of the flow of fermion number, in addition 
to the arrow on the external photon line associated with the A vertex. It is important to consider terms (f) and (g) 
separately, in order to show that the correct matrix element is produced by the pinching part of the kernel of the 
integral equation. The procedure is described in detail in Ref. [12 ] . 

The second diagram on the right hand side of Eqn. l[T4j) is of 1-loop order. Substituting the equation for D produces 
the graph in part (a) of Fig. [TH Substituting the equation for V produces the graphs in parts (d) and (e) of Fig. [TH 
Iterating the 1-loop term with itself produces the 2-loop diagram shown in part (b) of Fig. [T5l 




a bed 




FIG. 15: 2-loop contributions to A. 



Equation (l24l) ( as shown in Fig. [15]) is the same as the equation found in Ref. using the 2PI formalism (see 
Fig. 8 of Ref. [I4]), except that the vertices are now determined self-consistently through the integral equation ([TTll . 
The pair of equations (|24| (Fig. fTS"]) and ([TT]) (Fig. [4]) are precisely the two equations shown in Fig. 11 of Ref. [14 ]. 
and obtained previously, using kinetic theory, in Ref. [13 ] . 



VI. CONCLUSIONS 



We have shown that the 3-loop resummed 3PI effective action contains all of the leading order contributions to 
the conductivity. The two integral equations that resum pinch and colinear terms are produced naturally by the 
formalism, independently of any power counting analysis. The result is explicitly gauge invariant. The method we 
develop in this paper should be generalizable to the calculation of other transport coefficients, like the shear viscosity. 

Our calculation provides a connection between nPI effective theories and kinetic theories, and supports the use of 
nPI effective theories as a method to study the equilibration of quantum fields. 

This work was supported by the Natural Sciences and Engineering Research Council of Canada. 
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